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We present a semi-phenomenological treatment of light transmission through and its reflection 
from a ferrofluid, which we regard as a magnetically tunable system of dense random dielectric scat- 
terers with weak dissipation. Partial spatial ordering is introduced by the application of a transverse 
magnetic field that superimposes a periodic modulation on the dielectric randomess. This introduces 
Bragg scattering which effectively enhances the scattering due to disorder alone, and thus reduces 
the elastic mean free path towards Anderson localization. Our theoretical treatment, based on in- 
variant imbedding, gives a simultaneous decrease of transmission and reflection without change of 
incident linear polarisation as the spatial order is tuned magnetically to the Bragg condition, namely 
the light wave vector being equal to half the Bragg vector (Q). Our experimental observations are 
in qualitative agreement with these results. We have also given expressions for the transit (sojourn) 
time of light and for the light energy stored in the random medium under steady illumination. The 
ferrofluid thus provides an interesting physical realization of effectively a "Lossy Anderson-Bragg" 
(LAB) cavity with which to study the effect of the interplay of spatial disorder, partial order and 
weak dissipation on light transport. Given the current interest in propagation, optical limiting and 
storage of light in ferrofluids, the present work seems topical. 



I. INTRODUCTION 

Ferrofluids (magnetic nanoparticles dispersed in 
liquids) display very unusual magneto-optical properties 
(optical limiting, switching, and such like) that can be 
tuned by varying an externally applied magnetic field. 
In this work we report on first results of a combined 
experimental and theoretical study of light transmission 
through and its back-reflection from a magnetically 
tunable ferrofluid - a random suspension of magnetite 
(-^6304) nanoparticles in kerosene. This remarkable 
system permits spatial rearrangement of the particles 
magneto-statically so as to alter, via the structure factor, 
the effective scattering of light in the system. The dual 
role of the nanoparticles, acting both as dielectric 
Rayleigh scatterers and as permanent magnetic dipoles, 
is a salient feature of the nanofluidic optical system, 
inasmuch as it scatters dielectrically and may be ordered 
magneto-statically. 

The phenomenon of coherent multiple scattering of 
waves, like de Broglie (dB) electron waves or light waves 
in a disordered medium, has been extensively studied 
in the context of Anderson localization [U |2]. Here, 
in the limit of strong scattering, the elastic transport 



mean free path {1^) can decrease to the extent that the 
wave becomes spatially localized. This is the well-known 
loffe-Regel limit, fco-^e ~1, where ko (=27r/Ao) is the 
magnitude of the wave-vector corresponding to the 
wavelength Aq in the medium. In the case of electrons, 
strong scattering (£e < A^b) results in localization 
of the electron waves |3], as in the metal- insulator 
transition. For light, however, such strong scattering 
demands an unrealistically high dielectric (refractive 
index) mismatch, making it difficult to localize light 
[4j. At very long wavelengths, localization escapes 
the loffe-Regel condition because of the weakness of 
Rayleigh scattering (A,^'*). At very short wavelengths 
geometrical (ray) optics holds and there is no scope 
for localization, which is essentially a multiple wave- 
interference effect. The localization condition for light 
may, however, be closely approached by introducing a 
periodic modulation in the disordered medium such that 
the Bragg scattering opens a forbidden optical band-gap, 
presumably a pseudogap, giving a small £e ~ '^o due 
to enhanced random scattering. Indeed, this has been 
shown for photonic bandgap materials in the presence 
of disorder In our experiments, we combine disorder 
(dense random scattering) with partial order (Bragg 
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scattering) in the presence of weak dissipation obtaining 
in a ferrofluid placed in a magnetic field. We induce 
spatial order that can be magneto-statically tuned 
towards Bragg resonance so as to give rise to a much 
smaller ^e, causing the time spent by light within the 
medium to be significantly enhanced. In the presence of 
weak dissipation (inevitably present in a ferrofluid), we 
demonstrate, experimentally as well as theoretically, the 
simultaneous decrease of transmission and reflection of 
light over a range of disorder (£e-values). One may aptly 
refer to our ferrofluidic sample, essentially a weakly 
dissipative random suspension oi npF nanomagnets per 
unit volume, with a weak spatial periodicity superim- 
posed, as a Lossy Anderson-Bragg (LAB) cavity. Our 
dielectric nanomagnets scatter light by virtue of their 
complex (lossy) dielectric polarisability, presenting us 
with a model system that has in it an interplay of spatial 
disorder, partial order, and weak dissipation. 

In this work, we first present a brief discussion of the 
phenomenology of the long photon transit (sojourn) time 
that we expect in a ferrofluid due to multiple scattering. 
This is followed by ID numerical simulation of light 
transport based on invariant imbedding. We then 
present experimental results for our tunable ferrofluidic 
LAB cavity, which are found to be in general agreement 
with our theory. 



lengths in the LAB cavity than the transmitted ones; 
thus R diminishes to a smaller extent than T due to the 
dissipative effects. 

The above phenomenology remains essentially valid for 
a 3D system of dilute {^XqUpf < 1) random nanopar- 
ticulate (NP) Rayleigh scatterers (radius qnp <^ Aq). 
For such a system, the elastic and the inelastic mean 
free paths can be readily expressed in terms of the real 
and the imaginary parts of the optical dielectric constant 
{e = eji + iej ) of the dielectric nanospheres in the fer- 
rofluid. Explicitly, we have, = nppa, with the scat- 
tering cross-section, a, of the nanoparticles being given 
by 

1)2 

cr = —t^na-j^p- -^1:777 (2) 
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3 ''""^^(efl + 2)2 



and, = (2^e,)/(AoeJ^) 

In the case of our ferrofluid, however, we have the 
opposite limit of a dense random scattering system 
{^XqUpp > 1) with the number density npp of the scat- 
terers lO^^cm"'^ and wavelength of light. A, 633nm. In 
such a case, the scatterers lying within a wavelength scat- 
ter essentially in the forward direction, and coherently, 
calling for a different treatment as presented below for a 
ID model system. 



II. PHENOMENOLOGY 

Considering a LAB cavity of linear dimension ~ L, we 
examine the transmission through and reflection from it. 
The mean time tl for the photon to diffusively traverse 
the LAB cavity is given by 2Dtl=L^, with the diffusion 
constant, D — cqH^, cq being the speed of light in the 
medium. The actual path length traversed by the pho- 
ton is then A = Tpco. With inelastic scattering mean free 
path ii, the survival probability, Ps, for a photon injected 



into our LAB cavity becomes Ps 



It is this surviving photon that finally emerges from the 
sample - in transmission (T) or in reflection (i?) - for a 
ID system. For many multiple scatterings, an equipar- 
tition between the emergent quantities R and T is ex- 
pected, and consequently. 
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Clearly, R and T decrease as the ferrofluid is tuned to- 
wards stronger scattering. We note that, strictly speak- 
ing, Eq. 1 is valid for a system where the forward and 
the backward directions are clearly defined. Further, 
it pertains to the case where photons are injected well 
within the optical sample whose linear dimensions are 
much greater than £e- In a real experiment, of course, as 
also in the simulation to follow, photons are incident at 
one end of the sample. This leads to an asymmetry in 
that typically the reflected photons traverse smaller path 



III. ID SIMULATION BASED ON INVARIANT 
IMBEDDING 

We simulate our ID model system using invariant 
imbedding [6, 7 giving directly the emergent quantities 
R and T as functions of the sample length. Here we in- 
troduce a position (^)-dependent complex dielectric con- 
stant, e{t) that comprises three parts: a disordered part 
(efl(£)) which is random and real, accounting for elas- 
tic scattering; an ordered part, (ep(£)) which is periodic 
and real, representing Bragg scattering; and a dissipative 
imaginary part, iei{i) quantifying the loss. Also, eg rep- 
resents the real part of the medium's mean (background) 
dielectric constant, giving the wave velocity cq = c/y/eo, 
and the wavevector fco ~ ky/eo=2T:/ Xq. As is appropri- 
ate for a dielectric at optical frequencies, the magnetic 
permeability has been set to unity. Thus, we have for 
our ID model ferrofluid 



e{e) = eo + iei + eR.ie) + ep, 



(3) 



with 0< £ <L. We model efi{£) as a ^-correlated random 
Gaussian variable, and ep as a sinusoid of wavelength 
2tt/Q. The emergent quantities, R and T, can now be ob- 
tained as function of the sample length, L, by the method 
of invariant imbedding in the following equations that in- 
corporate all the essential features of our model: the real 
random disorder (?/p), the ordered periodic modulation 
(ryp), and the dissipation (777) : 

dR(x) i 

— ^ = 2iR{x) + - [1 + R{x f] [m{x) + Vpsin{qx) + ir]i] 
dx 2 

(4) 
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FIG. 1: (a) Plots of transmission jTp and reflection |_R|^ coefficients (using eqs. [4j [5| as function of the tuning parameter q 
for different values of the dissipation parameter (??/). Note the dips in both |r|^ and [Rf' at g = 2, corresponding to maximum 
scattering (Bragg condition). Here the disorder parameter t^q^ = 0.15, the depth of modulation J7p = 1, and the dimensionless 
sample length x = 7; (b) Plots of transmission |T|^ and reflection \R\^ coefficients (using eqs. Bl (sl as function of the tuning 
parameter q for different values of the disorder parameter {"Oor)- Note the optical stop band bracketting q — 2. Ifere the 
dissipation rji — 0, the depth of modulation rjp — 1, and the dimensionless sample length x = 14; (c) Plots of transmission 
\T\'^ and reflection \R\^ coefficients (using eqs. |4] [Hj) as function of the tuning parameter q for different depths of periodic 
modulation, rjp. Here the dissipation parameter rji = 0.2, the disorder parameter rjoR = 0.15 and the dimensionless sample 
length X — 7. The dominant role of the periodic modulation is clearly brought out here. 



and 

dTix) i 

— = iT{x) + -[l+R{x)T{x)][T]R(x)+ripsin{qx)+iT]j], 

ax 2 

(5) 

where the dimensionless variables and normalized param- 
eters are introduced as a; = kgL, q = Q/ko, r]ji{x) ~ 
efl(L)/eo, rip = ep/eo, and rji = e//eo- 

For the random variable r]ii{x), we have < rmix) >~0, 
< rifj(x)r]ji{x') >= TfQ^5(x — x'). Results of our numer- 
ical solution of the stochastic eqs. |4] and [5] are shown in 
Fig. 1 in the form of transmission (iTj^) and reflection 
(|i?| ) coefficients as function of the tuning parameter 
q for fixed sample length. This is shown in Fig. 1(a) 
for different values of the dissipation parameter 77/, for 
fixed disorder and fixed depth of periodic spatial mod- 
ulation. Fig. 1(b) shows the corresponding dependence 
on q for different values of disorder, for a fixed depth of 
modulation and zero dissipation, and Fig. 1(c) for differ- 
ent depths of periodic modulation. (A technical point to 
be noted here is that in our computation we have used 
the polar form for the complex quantities R and T, so 
as to simultaneously monitor both the magnitudes and 
the phases separately. This allows us to ensure the phys- 
ical upper bound of unity for both \R\^ and , for all 
realizations of disorder, dissipation and for all lengths). 
Note in Fig. 1(a) that when dissipation is not too small, 
both and \T\ dip as the Bragg condition [q ~ 2) is 
approached, with the dip being more pronounced for the 
former. Note also that \T\^ remains smaller than 



with the sum remaining less than unity (non-zero dissipa- 
tion). In Fig. 1(b), the opening of the pseudo-bandgap is 
clearly discernible (|i?|^ is nearly unity and |T|^ is nearly 
zero), with the sum remaining unity (zero dissipation). 
In Fig. 1(c) it is seen that periodic modulation plays a 
dominant role in the scattering. 

It is now instructive to compare these results of our ID 
numerical simulation with those from our experiments. 
Here we emphasize that the ID system is physically 
equivalent to a 3D system, with the proviso that all spa- 
tial fluctuations and modulations be taken along the di- 
rection of incidence of the light. This is indeed approx- 
imately the case for our ferrofluidic system inasmuch as 
the transverse magnetic field is expected to induce spa- 
tial ordering of the interacting nanomagnets, while the 
dense random dielectric scatterers scatter predominantly 
in the forward direction. Indeed, the Bragg scattering 
is the dominant determining feature as far as multiple 
scattering is concerned. 



IV. EXPERIMENTS 

We have measured the reflected and the transmitted 
intensities and their polarisations upon irradiation of our 
ferrofluid by a 633nm light from a 5mW cw He-Ne laser. 
The ferrofluid (Ferrolabs, USA) consisted of surfactant 
coated Fe^Oi nanospheres of radius ajvp (~ lOnm) ran- 
domly dispersed in kerosene. Figure [2] is the schematic 



4 



/ — \ 

D4 




FIG. 2: Schematic of the experimental setup, showing sample 
(S) kept between the pole pieces PI and P2 of an electro- 
magnet. G is the Gauss Probe, H the half-wave plate, D1-D8 
detectors, and Bl and B2 are non-polarising beam splitters 
while B3 and B4 are polarising beam splitters. Arrows B and 
k denote the directions of the applied magnetic field and the 
wave- vector of the incident light. 



of our apparatus. Linearly polarised light from the laser 
source (the plane of polarisation of which could be ro- 
tated by means of the halfwave plate, H) was incident on 
the sample (S) that was contained in a glass cuvette of 
size 1cm x \cm x 2mm (2mm along the direction of inci- 
dent hght). A non-polarising beam-splitter (Bl) placed 
in the path of the incidence allowed us to monitor the 
intensity of the incident light using detector Dl and that 
of the reflected light using detectors D5 and D6. The 
cuvette was held between the pole pieces (PI and P2) 
of an electromagnet that could produce a tunable mag- 
netic field tranverse to the direction of incidence of light, 
of upto 450G, uniform over the illuminated sample to 
within a few mG, as measured by the Hall-effect Gauss 
probe (G). The light transmitted through the sample was 
incident on a non-polarising beam-spiltter (B2); detector 
D2 provided a measure of the transmitted intensity while 
the combination of the polarising beam-splitter (B3) and 
detectors D3 and D4 enabled the determination of its 
plane of polarisation. Similarly, using polarising beam 
splitter B4 and detectors D5 and D6, the intensity and 
polarisation of the reflected light could be determined. 
Detectors Dl - D6 were identical photodetectors (Thor- 
labs DETllO). A small photodiode (D7) was placed in 
the narrow gap between the electromagnet and the cu- 
vette in order to measure the light intensity scattered to 
the side. A similar detector (D8) placed on top (per- 
pendicular to both the direction of light incident on the 
sample and to the applied magnetic field) was used to 
measure the light intensity scattered in that direction. 

Figure 3 shows the relative variation of the reflected 
and the transmitted intensities and their polarisations 
with applied magnetic field. In this case the incident 
light is horizontally polarised (i.e., plane of polarisation 
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FIG. 3: Oscilloscope traces showing, as function of time, the 
applied magnetic field and the light intensities as measured 
by detectors D2-D8. The inensity of the transmitted light re- 
duces to zero when the magnetic field is increased; that of the 
refiected light is found to decrease to a smaller extent. Both 
the transmitted and the reflected light retain the incident lin- 
ear polarisation, in this case horizontal polarisation. 

perpendicular to the incident wave-vector, k, and par- 
allel to the applied magnetic field, B). A simultaneous 
decrease in both the transmitted light (D2) and the 
reflected light (D5) is observed as the applied field is 
increased. While the transmission drops to near-zero, 
the reduction in the reflection is less pronounced. On 
decreasing the magnetic field back to zero, the original 
intensities are recovered. The intensities recorded on 
detectors D4 and D5 (horizontal polarisation) and D3 
and D6 (vertical polarisation) confirm that the transmit- 
ted and the reflected light retain the input (horizontal) 
polarisation. The scattered intensities in the transverse 
directions (D7, D8) are quite small. The situation 
was quite similar when the plane of polarisation of the 
incident light was perpendicular to the applied magnetic 
field (i.e., vertical). 

The experiment was then repeated for two other wave- 
lengths - (a) 612nm (He-Ne laser) for which the fall of 
the intensities with increasing magnetic field was slower 
than that for the 633nm wavelength, and (b) for 780nm 
(diode laser) for which too, the fall of intensity was 
slower than that for 633rtr7i light. As we shall see below, 
these experimental observations are all qualitatively 
consistent with what is expected from our ID model. 



V. DISCUSSION 

We now consider more explicitly how the applied mag- 
netic field tunes the scattering properties of the ferrofiuid. 
There is ample evidence that the nanoparticles (nano- 
magnets) rearrange to form spatially ordered structures, 
like linear chains directed along the magnetic field [MTT] . 
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An overall spatial ordering of the nanoparticles is to be 
expected on general grounds as resulting from the min- 
imization of the total energy of the magnetostatic cou- 
pling of the nanomagnetic particles with the externally 
applied magnetic field, and the magnetic dipolar interac- 
tion between these nanomagnets. The permanent mag- 
netic dipole moment of these single domain magnetite 
nanoparticles is 5 x 10** Bohr magneton. The magneto- 
static energy per nanoparticle in the external magnetic 
field of SOOG is about O.leV, exceeding both the thermal 
energy (ksT), and the inter-dipolar interaction magnetic 
energy per nanoparticle. This implies that the external 
magnetic field can indeed introduce spatial ordering of 
these nanomagnets. Thus, by varying the external mag- 
netic field, we can expect to continuously tune the spatial 
order (wave- vector Q), sweeping across the Bragg reso- 
nance condition {Q — 2fco, i.e., q = 2), leading to strong 
scattering. The strong scattering (small mean free path) 
enhances the time taken by light to transit through our 
LAB cavity. This should result in a dip in the transmis- 
sion and the reflection coefficients bracketing the Bragg 
condition, as indeed our simulations (Fig. 1) show and 
our experiments (Fig. 3) reaffirm. As noted before (fol- 
lowing Eq. [T]) the asymmetric decrease of the two in- 
tensities arises from the simple fact that light is incident 
from one end of the sample, and is not injected directly 
well within the sample (LAB). 

Encouraged by this qualitative agreement between the 
numerical results of our ID-model simulation and the ex- 
perimental observations on our 3D ferrofluid, we now pro- 
ceed to examine if the choice of the parameter values for 
the model (specifically, r]oR, rjp, rjj) do indeed correspond 
reasonably well to our experimental ferrofluid. The input 
parameters taken for the ferrofluid in the experiment are: 
nanoparticle mean number density npp — lO^^cm"'^, 
nanoparticle radius qnp — 15 nm, dielectric constant of 
the nanoparticle material Fe^Oi, e^p — 4-|-i2.2, and the 
dielectric constant of the suspension medium (kerosene) 
tK — 1-8. The above parameters straightforwardly give 
the mean dielectric constant (the real part used in Eq. |4] 
and|5] 

At: o n 47r o 

Co = £^(1 - -y'^ffO'np) + RcA<^Np)^aNP^FF - 2.1 

(6) 

Next, the dissipation parameter : 

47r o 

r]i = e//eo = Im{eNp)—nFFa]^p ~ 0.31. (7) 

The random disorder parameter tjor, arises out of the 
fluctuations of the local number density of nanoparticles. 
Assuming a Poissonian distribution for the fluctuations 
and making use of the fact that the variance equals the 
mean (both taken over the length scale of a wavelength) , 
we obtain 

rj'oR = 32n'iReep-e,n^a%pn^) 0.037 « 1, 

(8) 



as expected for a dense random system {^XqUff ^ 
1). This again is consistent with the fact that the dense 
random system is not, by itself, an effective scatterer - it 
is the spatial modulation (Bragg scattering) that makes 
the overall scattering effectively strong. 

The simultaneous minima of transmission and reflec- 
tion at g = 2 in Fig. 2 is a clear signature of the Bragg 
resonance condition, and is consistent with the ferrofluid 
parameters. This can be seen from the following 
geometrical considerations based on the fact that in a 
magnetic held (^ 100 G) the nanoparticles are known 
to re-arrange forming chains aligned parallel to the held 
[81-fTT|. Assuming that the overall number density of the 
ferrofluid remains unchanged, we have l/npF = 2aNpd?, 
where 2a mp is the nanoparticle spacing along the chains 
and d the interchain separation. This gives d ~half the 
wavelength of light, Aq in the medium, and thus closely 
corresponds to the Bragg condition. From the same 
geometrical considerations, the parameter rjp, giving the 
depth of the periodic dielectric modulation under the 
above Bragg condition, can be readily estimated. It is 
given by ryp = [I&tt /'i){eMP / eK){nFFa%ip) ~ 1. Overafl, 
given the complexity of the system and the simplicity of 
our model, we find the above estimates of the various 
parameters quite consistent with our experiments. 

Next we turn to polarisation, namely, how, despite 
multiple scattering, the transmitted as well as the re- 
fiected light retains the incident linear polarisation. Let 
us consider the case of transmission, where the initial 
and the final wave-vectors are parallel. For clarity, we 
resolve the linearly polarised light into its left and the 
right circularly polarised components, both of which tra- 
verse the same path-length in the optical medium, and, 
therefore, undergo the same dynamical phase shift. This 
causes no change in the state of polarisation. Phase 
shifts, however, do also arise as a result of the geometric 
phases suffered by the two circularly polarised compo- 
nents. These phase shifts are equal and opposite for the 
left and the right components, and have a magnitude 
given by the solid angle subtended at the origin in the 
space of directions traced out by the ray trajectory. For 
light, with photon spin =1, this does lead to a rotation of 
the plane of linear polarisation by half that solid angle, 
which is, of course, random for our disordered system. In 
the present case, however, the scattering is dominated by 
the Bragg scattering which is essentially in the forward 
and the backward directions. These multiple scatterings 
subtend zero solid angle, and hence give no rotation (due 
to the geometric phase) of the plane of incident linear 
polarisation for the transmitted (forward) and the re- 
fiected (backward) scattered light. Moreover, recall that 
a single scattering by a Rayleigh scatterer (our spherical 
nanoparticle) retains the incident linear polarisation of 
light, anyway. Besides, unhke the case in [TOl|Tl], our 
magnetic field is transverse to the direction of incidence, 
and hence there is no Faraday rotation. As already noted, 
in the transverse directions, the light intensities are very 
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small. 

Finally, our model can be viewed in a broader perspec- 
tive; some of the quantities of possible physical interest 
become calculable. It is possible to calculate the photon 
transit time (r/,) in our LAB cavity. Indeed, the intro- 
duction of a weak dissipation {£i) provides an effective 
measure |12) for tl, which can be readily shown to be 
given by e"^^^"/^' = + \R\^). Also, the amount of 
light energy stored {Us) in our LAB cavity under steady- 
state illumination (intensity /) can then be readily ex- 
pressed as Us = Itl- We note that this light is stored as 
light, and not as an electronic excitation. It is tempting 
to suggest here that this stored light may re-emerge from 
the LAB cavity as a flash if the magnetic field (and, there- 
fore, the confining strong scattering) could be switched 
off instantaneously. This, of course, assumes that the 
nanoparticles in the ferrofluid relax back sufficiently fast. 
The emergent light pulse would have a linear polarisation 
same as that of the incident light. Clearly, this would not 
be the case if the light had been stored as an electronic 
excitation, and re-emitted as fluorescence. 

VI. SUMMARY 

In summary, we have studied experimentally as well as 
theoretically light transmission through and its reflection 
from a magneto-statically tunable ferrofluid. We have 
developed a mechanism in terms of scattering of light in 
a disordered medium in the presence of weak dissipation, 
where the scattering could be effectively enhanced by a 



spatial periodic modulation of the background dielectric 
constant, giving Bragg scattering. Our theoretical 
treatment of transmission and reflection of light is based 
on the method of invariant imbedding. The key idea 
to emerge from our studies is that the strong multiple 
scattering of light and the resulting small mean free path 
can give rise to a long photon sojourn time, thereby effec- 
tively making our ferrofluid system akin to a lossy cavity. 
It is apt to point out here that the increased photon 
sojourn time calculated and reported in this work comes 
from multiple scattering because of disorder, re-inforced 
by Bragg scattering due to spatial modulation. It does 
not arise because of dispersion as in the well-known 
slow-wave systems, where the slowness comes essentially 
from the fast variation of the real part of the refractive 
index with frequency, as, for example, in the Elec- 
tromagnetically Induced Transparency (EIT) systems, 
where very low group velocities (as low as few cms~^ 
could be achieved over a narrow frequency window |13j . 
Given the recent research activity [101 IHl EH H] in 
the magneto-optical properties of magnetically tunable 
ferrofluids, the present work seems to be of some interest. 
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